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ON ! Abstract 

(2J[)' The systematic approach to study bound states in gluodynamics is presented. 

The method utilizes flow equations together with low-energy phenomenology, that 
provides the perturbative renormalization scaling in conjuction with the change of 
Cf~) ' the basis to constituent gluon states. The renormalized effective Hamiltonian of 

gluodynamics up to the second order is obtained at low energies, which provides a 
, kind of constituent gluon model for glueball bound states. The approach allows to 

J> \ include perturbative QCD corrections into nonperturbative calculations of many- 

body techniques. The performed numerical calculations support the constituent 
ly-^ . picture of hadronic observables. 

in : 
o . 

1 Introduction 

QCD is a widely accepted theory of strong interactions, where the phenomenon of asymp- 
totic freedom provides the success of perturbative calculations. There is still a gap be- 
tween perturbative behavior of QCD and its low-energy limit, where physical observables 
are described based on phenomenology. An ultimate goal of this study is to connect both 
regimes and obtain a low-energy effective Hamiltonian from canonical QCD Hamiltonian, 
^ ■ which then has to be solved for bound states. 

We suggest a framework which incorporates a perturbative behavior of QCD and our 
knowledge from QCD motivated phenomenological models. We consider pure gluody- 
namics in order to disentangle two problems of low-energy QCD; confinement and chiral 
symmetry breaking. 

Our primary aim is to obtain a renormalized gap equation in the gluon sector, which 
yields a gluon mass gap. The previous study concerns the fermion gap equation for dy- 
namical quark mass [0, 0. In this case a nonzero chiral condensate (qq) is produced, 
which leads to a spontaneous chiral symmetry breaking. In gluodynamics the massive 
gluon mode provides a double pole instead of a simple one in gluon propagator || . This 
is used as a phenomenological criteria of confinement |3|, f|, f|. Actually, we proceed 
in another way. We implement a confining potential, which provides in our calcula- 
tions nonzero gluon mass and gluon condensate (F^F^). One encounters a problem of 
UV-divergences with the Coulomb potential. For potentials which do not lead to UV- 
divergences, such as a pure confining potential 1/g 4 , unrenormalized gluon gap equation 
was obtained by Szczepaniak et.al. 0. First, we derive the renormalized gap equation, 
which provides a cut-off independent gluon dispersion relation up to leading order; it is 
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used then in the Tamm-Dancoff bound state equation to calculate glueball mass nonper- 
turbatively. 

We utilize the advantage of Coulomb gauge: the Coulomb propagator corrections give 
the complete QCD /3-function [0], @, |], that permits a simple implementation of renor- 
malization group improved perturbation theory. The Coulomb gauge allows an extension 
of the present approach to confining potentials as one uses in phenomenological analysis 

0- 

The motivation for this study is to set up a kind of constituent gluon model for 
glueball bound states, similar to constituent quark model. In the previous study, one 
made a Bogoliubov-Valatin transformation to a quasiparticle basis, with a BCS vacuum 
containing a qq condensate (pairing model) @. In the present work we apply the method 



of flow equations [T3] to transform to an effective Hamiltonian, which is block-diagonal 



in particle number space and describes quasiparticle excitations [[16 . 

The proposed approach utilizes the method of flow equations in theory of strong 
interactions with the confinement that is embedded. Using flow equations we scale down 
the (normal ordered) canonical QCD Hamiltonian from the bare cut-off scale A — > oo to 
some moderate scale A where the perturbation theory breaks down. Renormalization 
is performed by absorbing all (in this order of perturbation theory) UV-divergences 
in counterterm operators. This procedure recovers the perturbative renormalization for 



Hamiltonians [fl2|l , which aims to reduce the cut-off sensitivity of observables. Due to the 
principle of renormalization group invariance the physical gluon stays massless through 
the perturbative scaling. QCD coupling constant, renormalized to the third order, starts 
to grow, that stops asymtotic freedom at the scale A ~ A and forbids to use perturbation 
theory further. Strong coupling constant, a s ~ 1, does not reflect the whole complexity of 
strong interactions at low energies. We introduce confinement as a linear rising potential, 
that sets the (hadronic) scale for the gluon mass. In the renormalization group sense this 
"spoils" the theory: there arises the massive gluon mode, which breaks gauge invariance. 
In the present approach confinement makes it possible to run down flow equations in the 
"nonperturbative" region from Ao to the hadronic scale, where the effective Hamiltonian 
has block- diagonal form. The higher orders in iterrative procedure are suppressed by the 
inverse power of (heavy) gluon mass. By applying flow equations to block-diagonalize the 
"confined QCD" Hamiltonian one gets a closed chain of decoupled equations up to leading 
order, which can be solved analytically. The block-diagonal effective Hamiltonian, with 
fixed number of quasiparticles (constituent gluons and quarks), provides a constituent 
desription for hadronic observables |I8| . In the case of pure gluodynamics, one gets a 
constituent gluon model. 



2 Canonical Hamiltonian of gluodynamics in the Coulomb 

gauge 

The starting point is the canonical QCD Hamiltonian (pure gluodynamics) in the Coulomb 
gauge (V • A = 0) [|K]]. Physical degrees of freedom are the transverse gauge fields A and 



their conjugate transverse momenta II. The complete form of the canonical QCD Hamil- 
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tonian can be found also in M. There the potential source of difficulty comes from the 
Faddeev- Popov determinant J = det[V • D] (normalized by det[V 2 ] = 1), which arises 
when gauge fixing constraint equations for the physical fields are solved. We expand 
perturbatively the canonical QCD Hamiltonian in the bare coupling constant, keeping 
the terms to the second order. In leading order Faddeev-Popov determinant reduces to 
unity and simplifies the instantaneous term to have the pure Coulomb behavior. The 
canonical Hamiltonian in the gluon sector takes the form 

H C an = H + Hi n t • (1) 

The free gluon Hamiltonian, defined as H = H can (g = 0), is 

H = Tr| dx(n 2 (x)+Bi(x)) , (2) 

which is written through the abelian component of magnetic field B^, since the non- 
abelian magnetic field B = B a T a has components Bf = e ijk V ' jA a k + ^e ijk f abc A)A%. The 
trace in Eq. (Q) is understood in color space. The interaction part includes the (non- 
abelian) gluon terms: triple- and four-gluon vertices, and the Coulomb instantaneous 
interaction 

H int = Tr J rfx (B 2 (x) - B* (x)) + H c . (3) 

The Coulomb instantaneous interaction is given 

H C = -\ J rfxrfyp a (x)V(|x - y\)p a (y) , (4) 

with the Coulomb potential V(|x — y|) = — ot s /|x — y| (a s = g 2 /4Tr), and color charge 
density (only gluon component) is p a (x) = / afec A 6 (x)II c (x). 

We proceed in a standard way and express the Hamiltonian Eq. ([!]) in the Fock 
representation. The minimum ground state of the canonical Hamiltonian is achieved 
in the perturbative Fock space. Therefore, we choose the trivial (perturbative) vacuum 
|0) and construct the perturbative basis of free (current) particles: a^(k)|0) creates one 
(perturbative) gluon with zero mass, etc., and the vacuum is defined as a\0) = 0. We 
decompose the gauge fields 

= / A / i- [<(k) + af(-k)]e- kx 

J (271 ) 6 V2co> k 

n- (x) = -if ^/y[<(k) - «f(-k)K kx , (5) 

where the gluon energy = |k|. The canononical commutation relation in momentum 
space is 

[<(k), a^k')] = (2n) 3 S ab S^(k - k') £> y (k) , (6) 

where the operators a\ are represented in terms of polarization vectors as 
a"(k) = Yl\=x 2 e i(k, A)a a (k, A); and the polarization sum is 

Aj(k) = e i( k > ^) e i( k ' ^) = S ij ~ > ( 7 ) 

A=l,2 
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with the unit vector ki = ki/k. The condition of transversality is written 

k-a a (k) = k-a at (k) = 0, (8) 

and hence ki ■ Aj(k) = 0. 

We represent the Hamiltonian in the perturbative Fock space, Eq. @, and perform 
normal ordering in the trivial vacuum state |0). In addition to the canonical terms, normal 
ordering of the Coulomb interaction and the four-gluon vertex gives rise to one-body and 
condensate operators. 

The free gluon part Eq. (0) includes the kinetic term 



dk 



I- + u^af (k)<(k) + (— - ^4«(k)<(-k) + h.c.; 
a; k u; k 2 



f dk a+/1 1 r dk ,k 2 

J (kK(k) + 27 72^t-" k 



(2 

(af (k)<(k) 



(2vr) 3 V 

«(k)<(-k) + h.c.)) = H + (H - Ho) 



and the condensate term from normal-ordering of Eq. (|9|) 



Or 



dk ,k 2 



(2vr) 3 V 



+ w k 



where the volume is V = (27r) 3 <5^ 3 ^(0). In the perturbative basis, u; k 
Hamiltonian Eq. ([|) reduces to the canonical part 

dk 



Ho 



u^af (k)<(k) 



(2vr) ; 



(9) 

(10) 
the free 

(11) 



that defines the Fock space: Ho\n) = E n \n), where \n) is a set of eigenstates and the 
energy of the state \n) is given by a sum of one-particle energies E n = Y^m^m- 
The nonabelian gluon part includes in the order 0(g) the triple-gluon coupling 



l 9 jabc 



2^/2 



3 dk 

-j-J- I. IV;, 



k 



1.7 



(12) 



x: 



of(ki) + of f (-ki) I ^(kaJ+aJ'C-ka)! I<(k. 



where in the following ui\ = cj kl , etc. 

In the order 0(g 2 ) the normal-ordered four-gluon vertex 



H® 



x: 



a^+afi-^) a%k 2 ) + af(-k 2 ) 



i r <^k„ 



(27T) 3 ^(E k -) 



(13) 



af(k 3 ) + af( 



a*{k A )+af{-k A ) 



Normal ordering in Eq. ( |T3"D produces one-body operator 

dkd(\ 1 



n, 



(27r) 6 w k w q 



[2(J y - Ai(q)] 



(14) 



af(k)<(k) + 1 (<(k)a"(-k) + h.c. 
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and the condensate term 







a s 7r 



N C (N 2 C - 1)V 



dkdq 1 

(2tt) 6 w k a; q 



3 - (k)' 



(15) 



The Coulomb piece Eq. (H) includes the Coulomb instantaneous interaction 



(2) 



-fr^(n^)(2^< 3 >EM( 



^) V V(k 1+ k 2 ) (16) 



x : 



aj(ki) + af(-ki)] k(k 2 ) - of (-k 2 )l la, d (k 3 ) + af (-k 3 )l k e (k 4 



where V(k) = 47r/k 2 . The terms arising from normal ordering in Eq. (|16|) are 
the one-body operators 



n 



c 



+ 



4 

cy.<N r . 



dk dq 

(2tt) 6 
dk dq 
"(2^)6 



t/(k + q) 



V(k + q) 



£) y (q) af (k)a°(k) 



Aj(q) <(k)a. a (k) + h.c. 



(17) 



and the condensate term 



a 



a: 



c 



-^N C (N? - 1)V 



dk dq- 



7(k + q)P-l 1 + W • 



(18) 



The one-body and condensate terms diverge in the ultraviolet region. 



3 Hamiltonian renormalization 

As was mentioned in the previous section, normal-ordered canonical Hamiltonian includes 
UV-divergent terms. The first thing to do is to regularize these divergences (explicit form 
of regulator is specified below). The cut-off sensitivity is removed by renormalization. 
Renormalization is performed by adding counterterms, which are local operators with the 
symmetries of the canonical Hamiltonian. In order to find explicit form of the counterterm 
operators we use flow equations, which evolve the Hamiltonian from the bare cut-off scale 
A to some lower scale A A. It is enough to find the gradient of the Hamiltonian to 
define the counterterms. Also the form of the cutoff (regulating) function is specified by 
flow equations. 

In addition, new interactions, that do not change the number of particles, are gener- 
ated. The renormalized effective Hamiltonian is required to preserve the strucure of the 
canonical Hamiltonian. 

The renormalized Hamiltonian can be written in general 

H ien (A) = H can + 5X CT (A), (19) 
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where 5Xct{A) is a set of (unknown) counterterm operators. The renormalized Hamil- 
tonian depends implicitly on the cut-off A through the counterterms []. 

Primary aim of flow equations is to reduce the many-body problem of quantum field 
theory to a few-body one. The idea is to find a unitary transformation that transforms the 
Hamiltonian operator to a block-diagonal form, where each block conserves the number 
of particles. Generally, the Hamiltonian matrix can be represented in the particle number 
space as 

R _(PHP PHQ\ 

n ~ \ QHP QHQ J ' [ZL) 
where P and Q are projection operators on the subspaces with different particle number 



content. The flow equations [15] bring the Hamiltonian matrix Eq. (21) to the form 



_ / PH eS P \ 
HcS ~ { QH cS Q ) ' ^ 

where the two blocks of the effective Hamiltonian decouple from each other. It may be 
simpler then to solve for bound states within one block, say PH e ffP, than to diagonalize 
the complete Hamiltonian, Eq. (plD , of the original problem. Since, generally, the number 
of particles in P and Q spaces is arbitrary, one can reduce in this way the bound state 
problem with many particles to a few-body problem. 

The renormalization of possible ultraviolet divergences can be done also by the 
method of flow equations. By using flow equations to block-diagonalize the Hamilto- 
nian one eliminates the particle number changing contributions of the latter not in one 
step but rather continuous for the different energy differences in sequence. This procedure 
enables one to separate the ultraviolet divergent contributions and find the counterterms 
associated with these divergences. This covers the UV-renormalization for Hamiltonians 
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There is only one scale, introduced in the renormalized Hamiltonian by flow equations 
- the cutoff scale A, which defines the size of the canonical ("bare") Hamiltonian matrix in 
the "energy" space, and in the effective Hamiltonian plays the role of the regulator in the 



divergent loop integrals. This can be compared with the similarity scheme [12], where 
the effective Hamiltonian is band-diagonal in the "energy" space and one introduces 
therefore two scales - the A which is the size of the Hamiltonian matrix and the A C A 
which is the width of the band, and all the matrix elements of the effective Hamiltonian 
are squeezed in a band, i.e. \Ei — Ej\ < A. Generally, in the latter scheme not only the 
divergent contributions are regulated, but also any matrix element (its "external legs") 
is limited to permit the transition only between the states with the energy differences 
less than A. One should not associate any physics with the second (axilary) scale A, 

1 Actually, the A-dependence arises for the renormalized Hamiltonian, which is expressed in Fock space 
(A regulates UV divergent loop integrals in momentum space). In the given (perturbative, Wk = |k|) 
basis Eq. ( |l9| ) reads 

:ff rcn (A): = :H caa : + :SX CT (A): , (20) 
where ":" stands for normal-ordering in the (perturbative) vacuum. 
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and consider it as some technical tool to run the renormalization procedure and find an 
explicit form of the counterterm operators. 

The renormalized Hamiltonian, obtained by flow equations in the perturbative frame, 
preserves the form of the (original) canonical Hamiltonian. Therefore the renormalized 
Hamiltonian can be written in terms of field variables, regardless of basis, and can be 
decomposed in any nonperturbative Fock space, provided canonical form. Contrary to 
that, the effective Hamiltonian, obtained by the similarity scheme, can not be so trivially 
transformed to the arbitrary (nonperturbative) basis, since the effective Hamiltonian has 
the band-diagonal form only in the perturbative basis and does not preserve it elsewhere. 
This advantage of flow equations we exploit in section 4 to extend the calculations to 
the nonperturbative hadronic scales. 



The flow equations are written in differential form |y| 

dJ f- = w),m)\ 

r l (l) = [H i (l),H r (l)\. (23) 

where the Hamiltonian contains two pieces, H(l) = H^(l) + H r (l), with HJJ,) includ- 
ing particle- number conserving (diagonal) and H r (l) - particle-number changing (rest) 
terms of Hamiltonian; rj(l) is the generator of the unitary transformation Q; I is the flow 
parameter with the connection to the energy scale I = 1/A 2 , Iq = corresponds to the 
bare cut-off A — > oo. For the choice of the generator, given by Eq. fl25p, one ultimately 
eliminates the particle-number changing part of the Hamiltonian, i.e. H r (l — ► oo) = 
(actually the 'rest' part is exponentially suppressed for large I). In this approach the 
generator includes only the terms from particle number changing sectors. 

Decomposed in the perturbative basis (u^ = |k|) and normal ordered in the trivial 
vacuum |0), the canonical Hamiltonian, Eq. (||)-Eq. (|IBD, is regularized by the bare 
cutoff A — > oo, i.e. the Hamiltonian is written through the bare parameters (bare couling 
constants and masses) at Iq = 1/A 2 . The regulated canonical Hamiltonian is the starting 
point for the further calculations in section 3. 

Technically it is convenient to write flow equations for the coefficients before the 
operator structures. Flow equations are written selfconsistently for the coupling constants 
and masses (energies), that are the functions of the flow parameter I and, generally, the 
functions of the momenta: say, in Eq. ([5])-Eq. (|T8|) to the leading order the (triple-) 
coupling constant and the gluon energy become g(ki, k 2 , k 3 ); I and cj(k; I) as functions 
of in- and out- going momenta involved in the given sector. 

2 Unitary transformation reads 

HQ, l ) = U-\l, l = Q)H{l )U(l, l = 0) , (24) 

where 



t/(/,/ = 0)=Pexp U n( 



iVWj . (25) 
and P stands for ordering along the flow parameter /. The generator rj satisfies = —r\. 
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We evolve the Hamiltonian with flow equations to the second order. The energies get 
a correction in the order 0(g 2 ), i.e. cu(k; I) = |k| + 0(g 2 ). Therefore the free Hamiltonian, 



Eq. (|9p, is H = H + 0(g 2 ), with H given by Eq. flTT|). To the leading order we include 
triple-gluon vertex into the "rest" part H r and other terms of canonical Hamiltonian in 
the diagonal part Ha, i.e. 

H r = H^ + 0(g 2 ) 

H d = H + O(g 2 ) (26) 

where 0(g 2 ) stands for the second order terms. 

The triple-gluon vertex, obtained by symmetrization of Eq. (12), reads 



x: 



p (ki, k 2 , k 3 ; /)<(k 1 )a5(k 2 )a c fc (k 3 ) + 3^(k ls k 2 , k 3 ; I) a? (-k 1 )a$(k 2 )4(k 3 ) + h.c. 



with the form given by 

r ijfc (ki, k 2 , k 3 ) = - [(ki - h)j5 ik + (k 2 - fci)*Aj + (k 3 - k 2 )iS jk ] , (28) 
which satisfies 

r iife (-ki, -k 2 , -k 3 ) = (-l)ryjfe(ki, k 2 , k 3 ) 

T jik (k 2 , ki, k 3 ) = (-1)1^!, k 2 , k 3 ) , (29) 

In Eq. (p7|) go(l) and g\(l) are coupling constants defined further. To the leading order 
the generator of unitary transformation, defined as r]jp = [Ho, H^], is given 

tf> - ^/-/(n^)^)^(E^)^r yt ( kl , k2 , k3 ) 00) 



x: 



^(ki, k 2 , k 3 ; /)<(k 1 )a J 6 (k 2 )a c fc (k 3 ) + %(k ls k 2 , k 3 ; /)af (-k!) S 6 (k 2 )4(k 3 ) - h.c 



where 



?7o(ki,k 2 ,k 3 ;/) = D (^i, u; 2 , u; 3 )# (ki, k 2) k 3 ; I) 

r/i(ki,k 2 ,k 3 ;/) = D x {u u u 2 , u 3 )g x ()t u k 2) k 3 ; I) , (31) 



and the energy differences are 



D (ux,u 2 ,uj 3 ) = -(uji + to 2 + w 3 ) 

Di(a;i,a;2,a; 3 ) = -(-wi + cu 2 + u 3 ) , (32) 



In the order 0{g) the flow equation reads 



dl 
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[t#\£o], (33) 



that provides the solution for the coupling constants in Eq. ( p7| ) 



g (l) = g(0)e- D2 ° l 
g x {l) = g(0)e- D l l , (34) 

with (7(0) the bare coupling constant at A — > 00; we suppressed momentum idices in Eq. 
fl34|) . Equation ([34]) may be written in a general form [|17] 

&(0 = <KO)/(A;0> ( 35 ) 

where f(D;l) is the similarity function with the properties f(D;l = 0) = 1, f(D;l — > 
00) = 0. Such a more general form was first used by Glazek and Wilson []12[]. In Eq. 
fl34|) the gaussian similarity function is used; the sensitivity of bound state solution to 
the different choices of f(D;l) was considered in [l?]]. Corresponding to Eq. (|35|) , the 
generator of the transformation is written 

1 ( <Mf(Di;l) \ 

Index i shows, that coupling constant, <?j, and generator, 77^, are connected with the 
energy difference Di. We use this form, Eq. ( |35| ) and Eq. (|36|), in the calculations below. 

In the order 0(g 2 ) the flow equation for the particle number conserving part, H^, 
reads 

^ = [tfU^ (37) 

that contribute to the two-body, one-body and condensate (v.e.v.) operators in the 
Hamiltonian Q Further we consider all these sectors in sequence. 

3.1 Effective interaction (two-body sector) 

We calculate an effective gluon interaction in the color-singlet channel, that describes 
the glueball bound state a^a" 9 |0). For convenience we choose the cm. system P = 0. 
An effective interaction is written 

V eS = f^r'^J y^(k, q) :af (q)af (-q)a°(k)ag'(-k): ■ (39) 

where the interaction kernel V#/jy(k, q) includes the terms from the canonical Hamil- 
tonian and new interactions, generated to the second order by flow equations in the 

3 In the second order flow equation for the particle number changing terms, H r , reads 



iU =[r ] f\H^] + [r 1 f\H Q ], (38) 

(2) ~ (2) (2) 
where the second term with the generator ry g = [Hq, Hr ] makes H r to fall exponentially with flow 

parameter. 
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two-body sector Eq. fl3T|). Consider first the generated terms. The t-channel diagrams, 
arising from the commutator [a^a^a^, a^a^^a^y] (diagram without backward mo- 
tion) and from the commutator [a^a^a^, a^la^^a^,] (Z-graph in t-channel), do not 
contribute to the color singlet state; only s-channel diagram, coming from the commuta- 
tor [aiiCL2j a 3ki a i'i' a 2'i' a 3'fc']> nas t° be calculated (here the notation is = a^k^, etc.). 
In the cm. frame the flow equation for effective interaction is given 



dV gcn (l) _ ^jabcja'b'c f 

dl 8 J 



dkdq 1 



(27r) 6 ^ k ^ q ^ k _ q 

X T ijk (k, -q, -(k - q))r^/^(k, -q, -(k - q))D kk ,(k - q) ■ 4 • 9 

x [Vlfa, k, k - q; l)g x (k, q, k - q; I) + 7ft (k, q, k - q; l)gi(q, k, k - q; I)} 

x :af (^(-^(k)^; (-k): , (40) 

where factor '4' is a number of possible permutations; and property for T-factors Eq. 
( |29|) was used. The sum of two terms 'rjg' corresponds to the two different time-ordered 
diagrams in s-channel. We introduce the factor Sy^'ftf (k, q) as 

T ijk (k, -q, -(k - q))riY fc /(k, -q, -(k - q)) = g^Syj^y k , (k, q) . (41) 
which is given due to transversality condition Eq. (Bp 

Sijk,i'j'k'(k, q) = (kjdik - -{k + q)kSij + qi6jk)(kj>5i>k' - ^Ak + q)k'$Vj> + qi'Sj'k') , (42) 

this form is symmetric under transmutations k and q. The tensor structure of the effective 
interaction is defined by the contraction 
S ijkiiljlkl (k,q)D kk/ (k- q) = -%jy(k,q), that is given 

V ijt i>j>(k, q) = - yqiqi'Sjf + kjkj>6w + qi{kyb V j - kjb V y) + - kj>5ij) 



-(k — q) 2 ^niTti'5jjf + rtjUjidn' + mnjS^jf + UifTijfSij — niUjfSitj — Uiirijbi 



+ ^(l-W^I » ( 43 ) 

where rii = (k — q)j/|k — q| is a unit vector. Integrating flow equations Eq. (fl0|) with the 
generator given by Eq. ( |5B| ) provides the resulting generated interaction V gen (l — > oo) = 
V gen (the initial value is V^ en (^ = 0) = 0) 

abcf a> b > c fdkdq 1 T/ .fQ(D u D v ) , G{D 1 ,,D 1 )\ 



7 (2vr) b cu k u; q u; k _ q \ L>i / 

x lafCqJajVqXCkX/C-k):, (44) 
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where we have introduced the G-factor 

6( J D 1 , D 2 - 1) = - £ f(D 2 ; l')dl> 

&(D u D 2 ;l ^oo) = 9pi, A), (45) 
and the energy denominators are given 

D 1 = -(-a> q + w k + CU k -q) 

A' = -(-^k + + Uk-q) , (46) 

the interaction kernel Vy,i'j'(k, q) is defined in Eq. fl4"3|). For the 0-factors the following 
holds 

6(AA0 = ^(i-/ 2 (A0) 

e( A, A; + e(£> 2 , D i; = 1- /(A; 0/( A; 
e(D 1 ,D 2 ;l) , e(D 2 ,D r ,l) 1/1 1 



A ' A 2VA -xvvWC^O/CA;/)) 

1 ' ' 1 -V0( J D 1 , J D 2 ;/)-e( J D 2 , J D 1 ;/)), (47) 



2 VA A/ 

that we exploit further. The first equation exibits the meaning of B-factor: it plays the 
role of the ultraviolet regulating function in divergent loop integrals (see further). 

We considered the gluon interactions, generated in the second order by flow equa- 
tions, that conserve the number of particles. The generator r]jp Eq. ( |30D produces to 
the second order also the interactions, that change the number of particles, for example 
interaction W with the structure in second quantization o^a^a^a^. To eliminate these 
interactions we choose the generator of transformation r\ = r/^ + r] g 2 ', where the second 
order generator is calculated by ?^ 2 ' = [-ff , W 7 ] , that provides the exponential damping 
of the particle number changing interaction W. 

To maintain gauge invariance we add the normal-ordered Coulomb Eq. (|I~6"D and four- 
gluon Eq. ( |13D terms of the canonical Hamiltonian in the color-singlet channel. Projecting 
these terms on the color singlet state in cm. system gives 

vc = -f r bc r' b ' c J ^^^^ 

V 9 = ^ftofM'j ~ fyfrj) :af (q)a^(-q)<(k)<;(-k): , (48) 

Effective gluon interaction is given by a sum of terms 

V eS = V^ + Vc + Vg, (49) 



with generated interaction defined by Eq. ([441). 

Consider the leading behavior each of the terms Eq. fl49|) . It is convenient to introduce 

D = -(A + A'H-^k-q 

5D = lp 1 -D 1 ,) = -(w k -w q ), (50) 
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Actually the 0-factors in Eq. (fPj) are the functions of the dimensionless parameter 
£ = 5D/D |L7|1 . Therefore the combination in Eq. ((44[) is written 

— d— + — — = M 1 + -r^j- (51) 

where = 6(L>i, Dj,) - 0pi', £>i), and Q(D 1 , D v ) + Q(D 1 >,D 1 ) = 1. Generated and 
normal-ordered canonical interactions contribute to the effective interaction Eq. (|39|), 
respectively 



-%i'i'(k,q) 11+ 2 1 



2 y " v ^'^ v i-e ; 

«'#' = — F7v — 7a2^ + uj i) d ij d i'f 



(k-q) 

* 11' jy = y47r(5 lV %/ - , (52) 

For real processes (when energy is conserved, £ = 0) the interaction terms take the form 

ViTjj' = 17~ — %,i'j'(k, O) = if 7i V2 — ; 2?1 r\^»i»*'i'(^> ^) 

JJ 2 co£- q 2 (k — q) 2 + m 2 [k — q, A) 

W = "I (k-q)2 ^'' (53) 

with phenomenological dispersion relation cu 2 . = k 2 + m 2 (k, A) where m is dynamical 
gluon mass at the cut-off scale A, that we calculate further. The constituent gluon mass 
tends to zero as the scale is reduced and saturates to some value as the scale is increased 
Fig. (^). One can speculate about the leading order behavior of the effective interaction. 
The sum of Coulomb and generated terms behave as 

1 1 

VeS ~ J, ^ - J, V2~i J7\ TT , (54) 

(k — q) z (k — q) z + m z {k — q, A) 

that has 1/A; 2 singularity for the large energy A ^> lGeV and 1/k 4 singularity at low 
energy A < lGeV, corresponding to the Coulomb and the confining potentials, respec- 
tively. This can be considered as a selfconsistent condition to include confinement in the 
calculations below. 

3.2 Polarization operator 

Diagrams with two contractions of the commutator [rj^\ H^] contribute to one-body 
operators 

11 = ( n ?( k K 0t (k)a$(k) + <( k ) \ (<(k)4(-k) + h-c.)) , (55) 

where polarization operators n^(k) are given by the gluon loop integral 

^) = f^^M, (56) 
12 



and analogous for M^ b (k). Note, that while in general, from rotational invariance, one- 
body operators n?|(k) and M°f(k) have two components, 11$ (k) = A ab ( y k)5i j +B ab (k)k i k j 
(the same for M^ b ), only the term proportional to 8y survive in the integral Eq. ([j5|) 
due to transversality Eq. (^|). To separate A(k) one can choose the vector k = k z and 
consider components i = j = x,y. Then in Eq. (|56|), where n^k, q) has the terms of 

the type and q^j (there are no terms proportional to fcj and kj), one can substitute 

m - \{i - ikf) ■ (57) 

Considering only 5{j component in Eq. (|56|) simplifies the following calculations. 

Making use of the second order flow equations Eq. ([T|]), we have for the gluon polar- 
ization operator 



dIl ( l ) _ l N 5 aa> f dkdq 



dl 8 J (27r) 6 WkO'qCo'k-q 

X r ijfc (k, -q, -(k - q))r^ fc /(k, -q, -(k - q))D ii /(q)D fcfc /(k - q) • 2 • 9 • 2 

(/7i(k, q, k - q; /)#i(k, q, k - q; I) + r/ (k, q, k - q; l)g (k, q, k - q; I)) af(k)a"'(k) 
+ (»7i(k,q,k-q;0</o(k,q,k-q;0 + 7/o(k,q,k-q;Ofl , i(k,q,k-q;0) (58) 
x i«(k)<;(-k)+h.c.) 

where the first term (r]igi) corresponds to a gluon loop without backward motion, (r/ofi'o) 
stands for a Z-graph, and the last two terms (rjigo, Vo9i) are gluon loop diagrams with 
two incoming and two outgoing gluon lines; by calculating Eq. (]58f) we have also used 
the property Eq. (|29| ) for T-factors. 

The tensor structure of polarization operator is given by a double contraction of 
5(k, q) factor Eq. (|42|) with polarization sums in Eq. (^). We introduce therefore 
Fii#(k, q) = %,, v/ /{k. q)^ ) ji'(q)^ , fcfc'( k - q)> where 

F„(k,q) = 4,fcv((^ + ^ r ^ ¥ )(l-(^) 2 )) 

+ 2g ift , (l - 2(k P q)2 (l - (fcg) 2 )) , (59) 

Due to transversality condition we separate the 5a' component in polarization operator, 
Eq. (|5T|). The factor in Eq. ([59]) takes the form 



F«,(k,q)-><WG(k,q)/2 

2 4- „ 2 4- 

2(k-q) 2 



G(k, q) = 2(1 - {kqf) (k 2 + q 2 + n f Q \j l + (kq) 2 )) , (60) 



Integrating flow equation ( |58| ) gives 

<ik<iq 



5U 



H(0 - n(/ = 0) = a s7 riV c / G(k, q 
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' 9(A), A); &{Di,Di;iy 
. D Di . 



< + (k)<(k) 

' (61) 



where B-factor is defined in Eq. (fig) ; the energy denominators are given 

A = -(^ k + + ^k-q) 

A = -(-Wfc + ^ q + ^k-q) , (62) 

Quantum corrections have the opposite sign than obtained by the flow equations, since 
the flow equations carry out the scaling of Hamiltonian down from high to low energies. 
Therefore, the polarization operator, defined at some scale I = 1/A 2 , is given 

^ Oi s Ti „ T f dkdq 1 . 

n gen (A) = ^-N c / 7^—— G k,q) 

I J {ZTT) Wk^q^k-q 



X 



+ (i^ + ^) /(B ° ;A)/(Bi;A) + (^-i) (63) 



x (/ A ^^/(^; V)dY - f d ^f±f(D ; X)dX) ) i (<(k K (-k) + h.c.) 

we have used identities Eq. (H7|) to simplify this expression. 

In what follows we study the infrared behavior of the theory. In this case, external 
momenta are soft and A ~ A ~ — (^q + ^k-q) ~ — 2u; q in Eq. d63|). It is reasonable to 
assume for the regulating function / 2 (A A) = f(D;X), that corresponds to the rescaling 
of cut-off A/ v^2 — > A and does not change the physical result. The polarization operator, 
generated by the flow equations, takes the form 



2 



n gen (A) = -a s nN c f^-^^ 

J (27r) b CUk^q^k-q^q + ^k-q 

x (af (k)<(k) + ~«(k)<(-k) + h.c.)) . (64) 

We have used the free dispersion relation cj q = |q| in exponential factors f(D; A) up to 
the leading order. This is true for large cut-offs (see below). We add one-body operators, 
coming from normal-ordering of Coulomb and four-gluon vertex, Eq. (|D]) and Eq. (|14|). 
respectively. In these equations we also pick out the 5^ component using Eq. (|57D . The 
regulated polarization operators read 

n c (A) = ^Njp^-^V(k-0(l + (kr)e-^ 
8 J {2n) b co> k ^q 

x ((^ + ^) ?t(k)a«(k) + K-^)i(af(k)<(-k)+h.c.)) (65) 

2 J (27r) b Wk^q 
V(k)<(k) + i«(k)<(-k)+h.c. 
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We choose the regulating function in the loop integral over <iq in Eq. (p5|) in accordance 
with the regulator in the generated term Eq. ( |64D to match the energy denominator. 
Note that for the generated term the regulator arises naturally from the method of flow 
equations; for the normal-ordered canonical terms one needs to introduce the regulator 
in the loop integral in accordance with the generated terms. The resulting polarization 
operator is given 

n(A) = n gen (A) + n c (A) + n 9 (A) , (66) 



with the terms defined in Eq. (^) and Eq. fl6"5]). We use this expression in the next 
section to obtain gap equation. 

At the end of this subsection we calculate the counterterm, associated with the diver- 
gency of polarization operators. The generated, Coulomb and four-gluon terms contribute 
to the quadratically UV-divergent part, respectively 

™ = £^A'(-i + l + 2 

x / (^3^ (<'(kR(k) + i«(k)<(-k) + h.o.)) , (67) 
The corresponding counterterm is defined therefore 

SX f OT (A) = (a? (k)o?(k) + ~«(k)<(-k) + h.c.)) 

rfx:A a (x)A a (x):, (68) 



2 

m CT 



with the mass given 



m CT — ——Nc^A 2 , (69) 

7T 12 



When the quark sector is added using the same procedure, the algebraic coefficient in the 
propagator correction reproduces the QCD /3-function. This result was obtained in the 
similarity scheme [T2[] by D. Robertson et.al. [^]. This particular feature of the Coulomb 



gauge provides the simple implementation of the renormalization group improved per- 
turbation theory. 

Gluon condensate 

The second order flow equation Eq. ( |37j ) for the condensate term reads 



dO gen (l) l, r ,, r2 lVlJ . [ dkdq 1 



-AT c (iV c 2 - l)V / 



dl 8 C J (27r) 6 tUkCUqCJk-q 

xT ijk (k, -q, -(k - q))r^-/ fc ,(k, -q, -(k - q)) D iil (k)D jf (q)D kk/ ( 1 k. - q) • 6 • 2 
x?7o(k, q, k - q; l)g (k, q, k - q; I) , (70) 

where the only commutator [a^a^ja^, ai'l^f^k'] contribute to the vacuum expectation 
value; factor '6' is the number of possible permutations, and V is the volume. 
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The tensor structure of the condensate term is given by the triple contraction of the 
5-factor Eq. (f42|) with the polarization sums Sijk.i'j'k'Q*-, q)i}ji/(k)ZA,y(q)Z}fc£./(k — q) = 
G(k, q), where function G(k, q) is given in Eq. (|60|). Integrating flow equation ([T0|) 
produces 



O gcn (/) -O gcn (/ = 0) = 2a s7 riV c (iV c 2 - 1)V / ^ , ,\ ^Hj^o^j {n) 
The regulated condensate term is given 



dkdg 1 G(k,q) e(D ,D ; I) 
(27r) 6 ^ k ^ q u; k _ q 3 D 



O gen (A) = a 8 <KN c {Nt - 1)V / S ; /, G f n q) / 2 (^o; A) , (72) 

where Dq is given in Eq. (|62j), and because of the scaling down we changed the overall 
sign We summarize the regulated generated and normal-ordered, Eq. flI8|) and Eq. 
(|15D, condensate terms 



O gen (A) = -^N c {Nl-l)Y f^-^ ^ -( |q | + | k i + | k - q i)^ 

7 (27r) b 3a; k cc; q u; k _ q (cj k + cj q + ^ k _ q ) 

Oc(A) = ^/V (! (/V c 2 - 1)V / '^V(k - q) g - l) (1 + W) 2 ) e-'l-l+W) 1 ^ 
0„ (A) = _ 1)v / **_L. (3 - &)•) .-<^ , ,75) 

We have used the same reasoning as above to regularize the Coulomb and four-gluon 
terms. Radiative corrections to the gluon condensate, Eq. (|TUD, are summarized 

0(A) = O gcn (A) + O c {\) + O g (X) . (76) 

Calculate the counterterm associated with ultraviolet divergent condenstate part. The 
generated, Coulomb and four-gluon terms contribute, respectively 

div0 (A) = ^MzllX A . f-i + 1 + 2) / * J- . (77) 
V ; 7T 3 V 4 J J (27r) 3 2cj k V ; 



4 In the calculations below we use the symmetric form 



O gen (A) = a s KN c {Ni - 1)V 



(27r) 9 UJ1UJ2OJ3 



x (2n fS^ ( kl + k2 + ka) ^ f(D - A) , 



(73) 



where D is given in Eq. (|32|), and 

G(ki,k 2 ,k 3 ) = ([fc? + fc 2 2 + fc 3 2 ][l-2(fc 1 fc 2 )(fc 1 fc 3 )(fc2fc3)] + fc?(fcifc 2 ) 2 (fcifc3) 2 (74) 
+ fc 2 (fc 1 fc 2 ) 2 (fc 2 fc 3 ) 2 + fc 2 (fc 1 fc 3 ) 2 (fc 2 fc 3 ) 2N ) (2^) 3 5 (3) (k 1 + k 2 + ka) -» G(k,q) . 
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The counterterm is defined therefore 



SXUV = (^-l)v/^^ = ^/&(0|^(xK(x)|0> 1 (78) 

where the mass given in Eq. ( |69|) . Combining the counteterms Eq. (|68"D and Eq. (|78|), 

one gets 

SX CT = SX' CT + 8X£ T = m 2 CT Tr J rfxA 2 (x) , (79) 
which stands for the mass counterterm in the renormalized to the second order Hamilto- 



nian Eq. (19). Note, that there is one and the same coefficient ttict f° r both counterterms 
in one- and zero-body sectors. This suggests, that one can possible work in field theo- 
retical basis without decomposing the fields in Fock components (through creation and 
annihilation operators). 

At the end of this section we outline the Hamiltonian of gluodynamics, renormaized 
to the second order 

H r = H can + 5X C t = H + Rim + SX = Hq + R int , (80) 

where H can is given in Eq. ([!]) and the mass counterterm 5Xqt Eq. (|79| ) absorbes all 
possible UV-divergences, that appear to the second order. The second order perturbative 
corrections to this Hamiltonian include the condensate terms Eq. ([7^) and polarization 
operators Eq. (|66|). In one-body sector 



1 



(2tt) 3 k 1 v ' 1 w 2 J (2tt) 3V cu k 
(of (k)<(k) + i«(k)a«(-k) + h.c.)) , (81) 



2 

where m 2 CT is the mass counterterm Eq. (plf). In the 'leading log' approximation (for the 
terms propotional to A 2 ) the renormalized free Hamiltonian preserves the structure of 
the canonical free part Eq. (|TTD , provided the second order radiative corrections 0(g 2 ) 
are included, i.e. 

H r +O(g 2 ) = H , (82) 
that manifests renormalization group invariance (RGI), or coupling coherence in the 



context of Hamiltonian renormalization [ H[ . The same holds in zero-body sector for the 
condensate operator |. 



5 



The renormalized condensate term is given 



„ l /9 s,,/" e?k , k 2 + m^rp 



(2tt) s 

dk , ww k 



2 



(N? - 1)V / -^-c k + ^ 1)V C ' ct _ Wk) (83) 
where the second term vanishes (in the order A 2 ) when perturbative corrections are calculated. 



17 



One can proceed with flow equations order by order to find (all) the counterterms 
systematically. RGI insures, that the renormalized Hamiltonian preserves the form of the 
(original) canonical Hamiltonian, but only the coupling constants and the mass operators 
(that are usually classified as relevant and marginal operators in renormalization group 
sense) start to run with the cutoff scale. (We do not consider here, at least in the few 
lowest orders of PT, possible irrelevant operators, that may cause new type of divergences 
than are carried by coupling constants and masses). 

Using flow equations we run the renormalized Hamiltonian in perturbative basis, 
Eq. (^Of), downwards from the bare cutoff A to some intermediate scale A ~ ^qcd, 
where the perturbation theory breaks down. Due to the RGI the "physical gluon" stays 
massless through this perturbative scaling. We can not proceed with flow equations 
perturbatively further. The result of this stage is the renormalized to the second order, 
effective Hamiltonian, defined at some compositness scale, :H ren (A, A ): with A ~ Aqcd 
and A — > oo is UV bare cutoff, and semicolon means normal-ordering in the trivial 
vacuum |0). Though the renormalized Hamiltonian is obtained in the perturbative frame, 
it can be represented (regardless of the Fock basis) in terms of the fields A and II, Eq. 
( |88| ) . It is a consequence of the RGI. We denote the resulting renormalized Hamiltonian 
at the scale A as H ren (A,A ). 

4 Flow equations in the confining background 

We introduce confinement as a linear rising potential, that enables to run flow equations 
"nonperturbatively" until complete diagonalization of the Hamiltonian. In the renormal- 
ization group sense this " spoils" the theory: there arises the massive gluon mode, which 
breaks the gauge invariance. In the present approach the confinement makes possible to 
bring by flow equations the QCD Hamiltonian to a block-diagonal form with a fixed num- 
ber of quasiparticles in each sector. The elementary degrees of freedom (quasiparticles) 
become constituent gluons (quarks), which acquire masses of order GeV. Confinement 
(string tension) sets the (hadron) scale for the gluon mass. The effective (block-diagonal) 
Hamiltonian, written in terms of massive gluon modes, provides a kind of the constituent 
gluon picture. In the case of zero gluon masses one can not find 'sector' representation for 
the effective Hamiltonian, because of uncontrolable creation and annihilation of particles 
in vacuum. 

The instantaneous interaction, Eq. (f|), containes two pieces, the sum of the Coulomb 
and confining potentials 



in configuration space. The Casemir operator in adjoint representation is C a dj = (N% — 
1)/2N C , and a is string tension. Denote the renormalized effective Hamiltonian with the 
confinement embeded, Eq. (£53) and Eq. (TO), as H e ff(A,A ). 

As far as confinement is introduced the trivial vacuum |0) and the perturbative 
basis of free (current) particles, u>k = |k|, define no longer the "minimum" ground state. 
Therefore, we introduce the (arbitrary) basis, where the gluon energy is kept unknown, 



V(r) = -C a4r ^ + a\r 
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and is denned further variationally. Correspondingly, the (nontrivial) QCD vacuum \ 
is defined as a\Q) = 0, and the Fock space of constituent particles is given: a^O) creates 
the quasiparticle with the energy u^, etc. f] 

The renormalized effective Hamiltonian H e ff(A, Ao) at the scale Ao, written through 
the physical fields A and II and having confinement, is decomposed in the trial (con- 
stituent) basis and normal-ordered with respect to QCD vacuum |0). The resulting 
effective Hamiltonian is given by Eq. (|T2|)-Eq. (|T8| ) and Eq. fl81"|), Eq. (|83|), where the 
unknown gluon energy is variational parameter in the calculations. We denote the 
effective Hamiltonian in constituent basis as ::H e ff(A )::, where stands for normal- 
ordering in the QCD vacuum. 

We combine the terms in the effective Hamiltonian in each particle number sector 
according to the power of coupling constant 0(g n ) (n = 0, 1, 2). The higher order terms 
in the effective Hamiltonian are suppressed by the inverse powers of (heavy) gluon mass, 
which is of the order of hadron scale (see calculations below). Variational calculations 
below give for the gluon energy uj^ = |k| + m(k) with the effective mass m(k), where for 
small momenta m(k) tends to the value ~ GeV and vanishes at high momenta. Typically, 
in the second order of perturbation theory we have 

*W=V (87) 

where the matrix element of the canonical interaction V12 is of order of the inverse Bohr 
radius, ~ MeV, and the energy denominator, given by combinations of gluon energies, 
is ~ GeV. Though for small momenta the coupling constant is not small, a s ~ 1, we 
have the small parameter ~ 0.1 — 0.01, Eq. (^), due to the (heavy) effective gluon mass 
m(k = 0). 

In the absence of confinement the effective Hamiltonian preserves the form of canon- 
ical Hamiltonian due to RGI, with the change |k| — > properly. In the presence of 
confinement the canonical form is violated by the second order terms in the free Hamil- 
tonian, Eq. (PH), which contribute higher orders 0(g 3 ), etc. in flow equations. We aim 
to find the effective Hamiltonian after the scaling downwards from Ao to a hadron scale, 
say a/o". Since the effective Hamiltonian preserves the canonical form at least to the 
second order, the "perturbative" terms obtained by flow equations in section 3 match 



6 The change of basis from the (perturbative) current, Wk = |k|, to the (nonperturbative) constituent, 
with some w k , can be written as Bogoluibov-Valatin (BV) transformation from the "old", a, a 1 , to the 
"new", a,a\ operators 

a k = ch</> k a k + sh0 k aL k , (85) 



with BV angle 4>k given by ch(/> k = l/2(y / fc/oJ k + yu^Jk) The connection between the "old", |0), and 
the "new", \ fl), vacuum states is given 



1^) =ex p^E th ^ a k at - k ^ 10) 



(86) 

It was used in the work H to transform the QCD Hamiltonian into the constituent basis. 
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the "nonperturbative" terms arising when applied flow equations to H e ff(A, A ). Gener- 
ally, this approach allows to include perturbative QCD corrections into nonperturbative 
calculations of many-body techniques. 

We run flow equations in the confining background to block-diagonalize the effective 
Hamiltonian ::H e ff(A, Ao):: in the nonperturbative basis and to find consistently all the 
terms to the second order. Free Hamitlonian and the confining interaction are included 
in "diagonal" sector, the triple-gluon vertex forms " nondiagonal" sectors, that should 
be eliminated. The leading order generator of unitary transformation is given by Eq. 
(|3lD with the trial gluon energy u^. It is important, that the zero order approximation 
in flow equations include all particle number conserving terms of Hamiltonian. In the 
previous work |12[ free Hamiltonian was chosen to start iterative procedure. The better 



convergence to a bound state Hamiltonian is achieved in the case of flow equations. This 
enables to use an "effective" (trial) basis in calculations. At low energies this basis tend 
to constituent degrees of freedom. This description is close to the approach of Orsay 
group 0. 

We bring the Hamiltonian ::H e ff(A, A ):: to a block- diagonal form, where diagonal 
blocks decouple from each other including the second order. All calculations are the 
same as carried by renormalization procedure in section 3, except for the change of 
the perturbative basis to the unknown nonperturbative basis, fixed below (section 5) 
variationally. Since block-diagonalization is performed, the particle number changing 
interactions are eliminated (to the second order) completely. 

We outline the resulting block-diagonal effective Hamiltonian, renormalized to the 
second order, denoted as H e ff 

H e ff = H can + Vgen + 5X C T , (88) 

where H can coicide with canonical Hamiltonian Eq. ([!]) except for the particle number 
changing terms, which are eliminated in the leading order by the flow equations; in- 



stead the new gluon interaction V ge n is generated in two-body sector Eq. ([44]); the mass 
counterterm 5Xct Eq. fl79[) absorbes all possible UV-divergences, that appear to the 
second order. The second order perturbative corrections to this Hamiltonian include the 
condensate terms Eq. ( [75] ) and polarization operators Eq. 



Formally, the matrix elements of the effective Hamiltonian can be written 
(0|#eff|0) = O + O gcn (A) + O c (A) + O g (A) + 5X^ T 

(i|# eff |i) =h + n gen (A) + n c (A) + n 9 (A) + 5x' CT 

(2\H eS \2) = V gcn (X -> 0) + V C (X -> 0) + V 9 (X -> 0) , (89) 

where the term Oo, -f^o are defined in Eq. ( |T0"P and Eq. (||), respectively, the second 
order radiative corrections: condensate and polarization terms are in Eq. (|76|) and Eq. 
(ffijp, respectively, the countertems are in Eq. ([78] ) and Eq. QB8|) , the effective color- 
singlet interaction is in Eq. fl49|). Other possible two-body interactions, that change the 
particle number, are eliminated by flow equations to the second order and give nontrivial 
contributions to the higher orders. In the next section we use Eq. fl89|) to derive gap 
equation and Tamm-Dancoff bound state equation. 
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5 Glueball spectrum 



The result of the previous section is the renormalized effective Hamiltonian of gluody- 
namics up to the second order, Eq. Equations for physical observables based on the 
renormalized Hamiltonian are protected from UV-divergences. In this section we derive 
gap equation for the constituent gluon mass. Gluon dispersion relation and the effective 
gluon interaction are used to solve for the glueball bound state within the Tamm-Dancoff 
approach. 



5.1 Gap equation 

We use the traditional way to derive gap equation. In Bogoliubov-Valatin approach the 
trial vacuum state \Q) is optimized variationally, and the variational parameter is the 
Bogoliubov-Valatin angle of transformation. In our case the effective Hamiltonian has 
been transformed to a block-diagonal form in a quasiparticle basis. Hence the vacuum 
expectation value of the effective Hamiltonian should be optimized in this vacuum; the 
variational parameter is an unknown one-particle dispersion relation. One has therefore 

f W> = 0, (90) 
duj k 

that provides an equation for cj k . Making use of Eq. ( j89|) for the condensate terms the 
following gap equation is obtained 



2 i 2 i 2 

Lot = k + m 



CT 



- 2a,*Nj -X 1 fl + ^W q)e -V/A', (91) 



(27r) 3 W q ^k-q 

where D = — (cj k + co> q + co> k _ q ) and the mass counterterm m C T is defined in Eq. (|69|). 
To get Eq. (|9TD we have used for convenience the symmetrized form for the generated 
term, Eq. ([75]) . By performing variation over external gluon energy the regulator in the 
condensate terms reduces to the regulator, obtained in one-body sector for polarization 
terms when external momentum is soft, i.e. |k| |q|. 

The integral kernel of gap equation contains three terms. The most essential is the 
first one (with V), which forms the basis of the given above approach. We solve Eq. 
( pip numerically taking into account only the first term of the integral kernel and the 
corresponding to it Coulomb mass counterterm 

^t = --^A 2 - (92) 

Potential V is presented as a sum of Coulomb and confining potentials, Eq. (|90|), with 
the Fourier transform 

47T 1 

V(k)=C adj a s - + a- J (93) 
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where the latter interaction is added to obtain nonperturbative properties. The canonical 
Coulomb term (and also the neglected integral terms) leads to the leading A 2 and the 
subleading In A (corresponding to relevant and marginal operators in the context of 
renormalization group) UV-divergencies. These divergencies are regulated by the cut-off 
function R = exp(—q 2 /A 2 ). Since the gap equation is renormalized to the second order, 
the mass counterterm cancels the leading UV- divergency. Note, that confining potential 
also contribute a term of the type In A in UV region. The logarithmic cut-off dependence 
of the Coulomb and the confining integral terms leads to a slow logarithmic growth of the 
mass gap with A. This behavior can be absorbed by the running of the strong coupling 
constant a s . 

Though the integral kernel of gap equation is obtained to the second order in per- 
turbation theory, the numerical solution of this equation provides the gluon mass non- 
perturbatively. The linearized gap equation is solved for u;(k). The dispersion relation 
is displaced in Fig. (]]]). The free behavior cj(k) = k is recovered at high energy, and at 
low energy uj tends to the constituent gluon mass of roughly 0.9 GeV. This suggests the 
suitable parametrization cj(k) = k + m(k), which defines the 'running' gluon mass m(k), 
Fig. (|2|). A good fit to the numerical solution of the linearized gap equation is obtained 
by 

k 

m(k) = m(0)exp( — ) , (94) 

K 

with an effective mass m(0) = 0.9 GeV and k = 0.95 GeV. 

The gluon condensate is another nonperturbative characteristic calculated within this 
approach. The condensate is given by the vacuum expectation value (0|n 2 + B 2 |0) with 
only abelian component for the magnetic field (pairing ansatz), Eq. (0). We regulate 
this value by subtracting the perturbative contribution 

a s NZ-1 r°° 2 Mk) - kf 

{ V F ^ } -—^i dkkas Mk) ' (95) 

Using the dispersion relation obtained above ui(k) the gluon condensate is obtained 
1.3 • 10 _2 Ge^ 4 (for the cut-off A = 4 GeV), that agrees with the sum rules [22]. The 



dependence of the condensate on the cut-off is shown in Fig. (|J). The behavior of the 
effective gluon mass m(0) with A is displaced in Fig. (||). Logarithmic dependence of 
both terms can be eliminated by including the renormalization group running of coupling 
constant. 

5.2 Tamm-Dancoff approach for glueball spectrum 

We approximate a glueball bound state to consist of two constituent gluons. The glueball 
wave function in the rest frame is 

IVO = / ^^'(q)af(q)«f (-q)|o> , (96) 
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where is a creation operator of quasiparticle with an effective dispersion relation u; q . 
Tamm-Dancoff approximation works reasonable for a block- diagonal effective Hamilto- 
nian. To obtain the Tamm-Dancoff bound state equation we project the Schrodinger 
equation H\ip n ) = E n \ip n ) on the two-body sector; the result reads 



(E n — Eq 



*&(q) 



(q, 



with the notation 



q)) = Ai'(q)A-i'(q)^"(q) + A;'(q)Av(q 



-q 



(97) 



(98) 



and ajj 



a*+(q) 



color indices are omitted. In Eq. fl98f) we subtract the trivial vacuum 
energy i?o ; defined as i?|0) = Eq\0). Tamm-Dancoff equation with all possible terms 
obtained to the second order is given in Appendix (|A]). We consider here the most 
essential part of the integral kernel, neglecting the perturbative transverse gluon exchange 
and terms coming from the normal-ordering of four-gluon vertex. The contribution of the 
neglected terms to the mass spectrum is expected to be small. Tamm-Dancoff equation 
for (pseudo) scalar glueball states (with total angular momentum J, parity P, and charge 
conjugation C) reads 



(E n - E )<j) n (q) 
dp 



'q 2 + m 2 CT 



a 



— N, 



(2tt) 3 
dp 



K(p-q) 
V(p - q) 



-p 2 /A 2 



UJpUJq 



4>n(q) 



2(pq)F JPC (p,q)Mp) 



UJ p UJ n 



(99) 



with 



ip>qj 
(p>q) 



i + 

2 



(100) 



and the Coulomb counterterm m given by Eq. fl92|). 

Bound state equation ( ^8|) has two types of divergences, UV associated with canonical 
Coulomb interaction and IR coming from confinement. UV divergences in the self-energy 
term are regulated by the cut-off function, and the leading divergent part is canceled by 
the mass counterterm m. UV behavior of potential part (the last term) is regulated 
by the wave function. Though kinetic (self energy term) and potential parts contribute 
the infrared divergent pieces the bound state equation is infrared finite due to complete 
cancelation of both parts in infrared region. This cancelation happens for the color-singlet 
objects 0. 



Numerical calculations of Tamm-Dancoff equation (|99| ) are performed variationally 
with a set of gaussian test functions. Results of calculations for the lowest glueball states 
are presented in Table ([!]) and compared with the available lattice data [p3[ . 

Lattice calculations are done for SU (3) pure gluodynamics, using anisotropic lattice 
and improved SII action. Better agreement with the lattice data is achieved for the scalar 
channel. Remarkable, the mass of the lowest scalar glueball ++ is roughly twice of the 
effective gluon mass m(0) obtained before. 
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0++ 


*++ 


o-+ 


*-+ 


Tamm-Dancoff, (MeV) 


1760 


2697 


2142 


2895 


lattice data J23fl, (MeV) 


1730(80) 


2670(130) 


2590(130) 


3640(180) 



Table 1: Glueball spectrum for the lowest scalar and pseudoscalar states (a s = 0.4, a = 
0.18GeV 2 , A = AGeV, N c = 3). 

6 Conclusion 

We outlined the systematic procedure to solve the bound state problem in QCD. Our 
study is driven by the success of the constituent quark model. The approach is based on 
the method of flow equations, that includes the renormalization group scaling of canonical 
Hamiltonian from high to low energies and simultaneous change of basis from current to 
constituent degrees of freedom. The renormalization group aspect of the method protects 
the observables, calculated from effective Hamiltonian to be practically free from UV 
divergences. Therefore we are able to include radiative corrections at one-loop level. 
The leading UV behavior is canceled by the mass counterterm. The slow logarithmic 
dependence of masses on the cut-off, which is still left after perturbative renormalization, 
can be absorbed by the nonperturbative renormalization of strong coupling constant. 
Confinement interaction can cause the infrared problem. In glueball bound state equation 
this potential problem is avoided by adding a diagram with soft gluon emission and 
absorbtion in the second order. This leads to complete cancelation of infrared divergent 
parts for color singlet objects. 

The use of effective basis (built by including all possible particle number conserving 
interactions in Hd and acting in QCD vacuum) in conjunction with the renormaliza- 
tion group scaling makes possible to get an effective block- diagonal Hamiltonian with 
fixed number of quasiparticles (constituent gluons). An effective Hamiltonian provides a 
constituent description for glueball bound states at low energies. 

Based on flow equations, we are able to extend the study to nonperturbative range of 
energies. The renormalized effective Hamiltonian, obtained in perturbative frame by flow 
equations and including confining interaction, can be represented in the nonperturbative 
(effective) Fock space provided the canonical form. Contrary, the effective Hamiltonian, 
obtained in similarity scheme, is band-diagonal only in perturbative frame and can not 
be transformed in nonperturbative basis. Though coupling constant is not small at low 
energies, a s ~ 1, we are still able to perform analytical calculations by flow equations, 
since higher orders are suppressed by a (heavy) effective gluon mass. 

To provide a detailed analyses and to compare with lattice data we solve numerically 
the gap equation for an effective gluon energy and Tamm-Dancoff glueball bound state 
equation. Roughly the glueball mass is twice of the effcetive gluon mass, that supports 
the constituent picture. The gluon condensate based on the effective gluon dispersion 
relation is in agreement with QCD sum rules. 
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q(GeV) 

Figure 1: One particle dispersion relation. Dots represent the numerical solution of gap 
equation u;(k) (a s = OA, a = 0.18GeV 2 , A = AGeV, N c = 3), the solid line stays for the 
free dispersion relation cu(k) = k. 
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Figure 2: Gluon mass. Dots represent the numerical solution for m(k) = u(k) — k (a s = 
0.4, a = 0.18GeV 2 , A = AGeV, N c = 3), the solid line is a fit m(k) = 0.9 * exp(-A;/0.95) 
(parameters are in GeV). 
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Figure 3: Cut-off dependence of the effective gluon mass (a s = 0.4, a = 
0.18GeV 2 , N c = 3). 
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Figure 4: Gluon condensate (a s = 0.4, a = 0.18GeV 2 , N c = 3). 
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A The complete form of Tamm-Dancoff equation 

We outline the complete form of Tamm-Dancoff equation with all possible terms obtained 
to the second order. One- (only a^a component) and two-body sectors, specified in Eq. 
(j8~9D , contribute to the color singlet glueball state. Making use of Eq. (|97|), we get the 
resulting Tamm-Dancoff equation for glueball bound state 



(E n - Eo) (Xg(q)^(q) + ^'(q)^(q 
dp 



'q 2 + m 2 CT 



+ Wo 



+ a s 7rN c 



+ a s irN c 



(2tt) 3 
dp 



Hp <l) (!-' 
1 



(27r) 3 ^ p ^ q 

dp 1 



(¥fK;A) 



G(p,q) 



(2tt) 3 w p ^ q u; p _ q 

x te(q)0f(q)+^'(q)0-(q) 



7 2 (A;A) + / 2 (A;AV 



Do 



A 



-iV r 



- ct s 27riV c 



(2tt) 3 
dp 2 
(27r) 3 a;pa; q 
dp 



(jf2»(q)xS»(p)0t(p) + ^(q)^ m (p)0- n (p 



Ai(q)An( P ) - i ^(q)Xr( P ) 
/0(A,A') , 6(A',A)\ 



(27r) 3 WpW q w p _ q \ A A' y 

(q, P ) (xg,(q)x t r'(p)^(p) +^(q)^r'(p) 

where the energy denominators are 



tin / 



(P) 



(101) 



-Do = -(W q +Wp + W p _ q ) 
A = -(-OJ q + Up + U; p _ q ) 
A' = -(-W p + ^ q + Wp-q) , 



(102) 



and the following notations were introduced 



tin / 



IP) 



dn i 



1 

2 1 

1 



^(p)+0 /n (- P )) 



dn i 



tin/ 



-P)) 



*£(q) = Ai(q)^n(q) + An(q)^i(q) 
y£(q) = Ai(q)An(q) " An(q)%(q) 



(103) 



Canonical UV divergencies are regulated by the cut-off function, that is specified in the 
main text; factor Sk m> k'm' is given in Eq. (|42|), and 9 factors in Eq. (|45|) . We write Tamm- 
Dancoff equation explicitly in the scalar (0 ++ ) and pseudoscalar (0~ + ) channels. For ++ 



28 



the wave function can be written </> ij = 5 l ^(j){q). Tamm-Dancoff equation takes the form 

'q 2 + m 2 CT 



ljj r 



+ o; q 



+ a s irN, 

+ a s nN, 
+ 

+ a s 7rN c 



dp 

(2^ 
dp 1 

(27r) 3 o;pa; q 
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(2tt) 3 cj p u; q u;p_ q 



3 ./> P ;A) 



G(p,q) 



7 a (A>;A) , / 2 (^i;A)' 



a s 2nN c 



(27T) 

<ip 1 
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<ip 1 



(3 - (m) 2 )^(p) 



(27r) 3 a;pO; q a;p_ q \ L>i 



e{D Xi D y ) Q{D V ,D X ) 



G(q,p)0 n (p) 



, (104) 



In the pseudoscalar channel, + glueball state, the wave function can be represented 
u (q) = e^ k q k (f)(q). Tamm-Dancoff equation reads 



(E n - E )4> n (q) 
dp 



'q 2 + m 2 CT 



+ a s irN c 
+ a s 7cN c 



(27T) 
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WpWq 



2(pg)0 n (p) 



(27r) 3 a;pa; q cc;p_ q \ L>i 



x 2 P V(i-(pqr) - 



i 



vq [p - q) 



(pq) 0n(p) 



(105) 



The reduced version of Eq. ( |104| ) and Eq. ( |105| ) with the Coulomb term in the integral 
kernel is used in the main text for numerical calculations. 
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